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ABSTRACT 

We present calculations of r-process nucleosynthesis in neutrino-driven winds from 
the nascent neutron stars of core-collapse supernovae. A full dynamical reaction network 
for both the a-rich freezeout and the subsequent r-process is employed. The physical 
properties of the neutrino-heated ejecta are deduced from a general relativistic model in 
which spherical symmetry and steady flow are assumed. Our results suggest that proto- 
neutron stars with a large compaction ratio provide the most robust physical conditions 
for the r-process. The third peak of the r-process is well reproduced in the winds from 
these "compact" proto-neutron stars even for a moderate entropy, ~ 100 — 200A^a^, and 
a neutrino luminosity as high as ~ 10^^ ergs s~^. This is due to the short dynamical 
timescale of material in the wind. As a result, the overproduction of nuclei with A <^ 120 
is diminished (although some overproduction of nuclei with A 90 is still evident). The 
abundances of the r-process elements per event is significantly higher than in previous 
studies. The total-integrated nucleosynthesis yields are in good agreement with the 
solar r-process abundance pattern. Our results have confirmed that the neutrino-driven 
wind scenario is still a promising site in which to form the solar r-process abundances. 
However, our best results seem to imply both a rather soft neutron-star equation of 
state and a massive proto-neutron star which is difficult to achieve with standard core- 
collapse models. We propose that the most favorable conditions perhaps require that a 
massive supernova progenitor forms a massive proto-neutron star by accretion after a 
failed initial neutrino burst. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: abundances 
— stars: mass loss — stars: neutron — supernovae: general 
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1. INTRODUCTION 

The r-process accounts for the origin of about a 
half of the abundances of elements heavier than iron. 
The other half mostly come from the s-process. Nev- 
ertheless, the astrophysical site for the r-process has 
not yet been unambiguously identified. Among the 
large number of proposed candidates, the neutrino- 
heated ejecta from a nascent neutron star ( hereafter 
"neutrino-driven wind" ) has been sug gested ( Woosley 
k Hoffman 1992 ; Meyer et al. 1992 ) as perhaps the 
most promising site. Woosley et al. (1994) demon- 
strated that the solar r-process abundances were well 
reproduced as material was ablated from the proto- 
neutron star in neutrino-driven winds. There are, 
however, a few serious problems in their numerical 
results. First, elements with A ^ 90 were signifi- 
cantly overproduced by over a factor of 100. Sec- 
ond, the requisite high entropy (^ 400A'Afc) in their 
supernova simulations has not been duplicated by 
other independent theoretical studies (Witti, Janka, 
& Takahashi 1994; Takahashi, Witti, & Janka 1994; 
Qian & Woosley 1996). On the other hand, other 
viable sites have been demonstrated to also natu- 
rally reproduce the solar r-process abundance pat- 
tern. For example, it has been shown that neutron- 
star mergers can naturally reproduce the abundances 
of nuclei with A > 130 (Freiburghaus, Rosswog, & 
Thielemann 1999). Collapsing 0-Ne-Mg cores result- 
ing from progenitor stars of ^ 8 — IOMq has also been 
suggested (Wheeler, Cowan, & Hillebrandt 1998) to 
be a promising site for the r-process. 

In contrast to the above difficulties in the theo- 
retical studies of the neutrino-driven wind scenario, 
observational data seem to confirm that Type II su- 
pernovae are indeed the site for the production of r- 
process elements. For example, spectroscopic studies 
of metal-poor halo stars indicate very early enrich- 
ment of r-process elements consistent with production 
in massive stars (Sneden et al. 1996, 1998; Ryan, Nor- 
ris, & Beers 1996). Using a chemical evolution model 
Ishimaru & Wanajo (1999) have furthermore demon- 
strated that the large dispersion of the r-process ele- 
ment europium in metal-poor halo stars is also repro- 
duced if the r-process originates from Type II super- 
novae (see also Tsujimoto, Shigeyama, & Yoshii 2000; 
Qian 2000). 

As one way to fix the schematic wind models, 
it has been suggested that general relativistic ef- 
fects increase the entropy and reduce the dynami- 



cal timescale of neutrino-driven wind models. Us- 
ing semi-analytic studies for spherically symmetric, 
steady flow of neutrino-heated ejecta, Qian & Woosley 

(1996) showed that the inclusion of a first post- 
Newtonian correction to the gravitational force equa- 
tion increased the entropy by ^ 60% and reduced 
the timescale by a factor of '--^ 2. Cardall & Fuller 

(1997) further developed this argument by consider- 
ing a fully general relativistic treatment. They consid- 
ered a wide range of neutron-star compaction ratios 
(ratio of the gravitational mass to the neutron star ra- 
dius in Schwarzschild coordinates). They showed that 
a more compact neutron star leads to significantly 
higher entropy and a shorter dynamical timescale in 
the wind. In their study, however, the sensitivity to 
the neutron star mass and neutrino luminosity were 
not considered. 

Otsuki et al. (2000, hereafter Paper I) studied 
physical conditions of neutrino winds in a manner 
similar to that of Cardall & Fuller (1997), for wide 
ranges of the neutron star masses and the neutrino 
luminosities. They suggested that high entropy and 
short dynamical timescales were obtained for a mas- 
sive neutron star and a high neutrino luminosity. 
However, the radii of the neutron stars were fixed 
to be 10 km in their study. Moreover, these stud- 
ies did not perform calculations of realistic r-process 
nucleosynthesis except for one specific case of a neu- 
tron star mass (= 2Mq) and a neutrino luminosity 
(Li,_j — 10^^ ergs s~^) in Paper I. 

The purpose of the present study is, therefore, to 
examine more quantitatively the r-process in neutrino- 
driven winds utilizing a fully implicit reaction net- 
work of over 3000 isotopes. Trajectories of thermo- 
dynamic quantities for material in the winds are de- 
rived for various combinations of three parameters: 
the neutron-star mass; radius; and neutrino luminos- 
ity. The semi-analytic, general relativistic model de- 
veloped in Paper I is adopted to describe material in 
the neutrino-driven wind. 

The paper is organized as follows: In the next sec- 
tion, the dependences of the entropy and the dynam- 
ical timescale on the neutron-star mass, compaction 
ratio, and neutrino luminosity are discussed. In § 3, 
the r-process nucleosynthesis is performed using the 
neutrino-driven wind trajectories obtained in § 2. 
The yields for various neutrino luminosities are mass- 
integrated by assuming an exponential time evolution 
for the neutrino luminosity. These are compared to 
the solar r-process abundances. The implications of 
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Table 1 
Model Parameters 





A 


B 


C 


M/Mq 
R (km) 
GM/c^R 


1.4 
10 
0.21 


1.4 
7 

0.30 


2.0 
10 
0.30 



this study are discussed in § 4. 

2. NEUTRINO-DRIVEN WINDS 

The models of neutrino-driven winds used in this 
work were developed in detail in Paper I. We briefly 
describe the models here and point out some improve- 
ments which have been added. The system is treated 
as time stationary and spherically symmetric. The 
physical variables in the neutrino wind are then func- 
tions of the radius r only. The ejected mass Mcj 
by neutrino heating is assumed to be negligible com- 
pared to the mass of the neutron star M . Therefore, 
the gravitational filed in which the neutrino-heated 
matter moves can be treated as a fixed background 
Schwarzschild metric. The equations of baryon, mo- 
mentum, and mass-energy conservation then become. 



M 
du 
dr 



A'Kr pu , 
1 + {u/cf 



2GM/c^r dP GM 



de 
dr 



P dp 
p^ dr 



P/c2 dr 



(1) 
.(2) 

(3) 



where M is the mass ejection rate from the surface of 
the neutron star, q is the heating rate, p is the mass 
density, P is the pressure, and e is the specific internal 
energy (Shapiro & Teukolsky 1983). The velocity u 



is related to the proper velocity of the matter v as 
measured by a local, stationary observer by 



[l + iu/cf ~2GM/c^r] ^''^ u 



(4) 



The source term q includes both heating and cool- 
ing by neutrino interactions. Heating is due to neu- 
trino and anti- neutrino captures by free nucleons, 
neutrino scattering on electrons and positrons, and 
neutrino-antineutrino pair annihilation into electron- 
positron pairs. The redshift of neutrino energies and 



bending of the neutrino trajectories by general rela- 
tivity are explicitly taken into account. Cooling is due 
to electron and positron capture by free nucleons, and 
annihilation of electron-positron pairs into neutrino- 
antineutrino pairs. For the latter, the more accurate 
table of Itoh et al. (1996) is utilized rather than equa- 
tion (15) in Paper I. All other neutrino heating and 
cooling rates are taken from Paper I. 

The neutrino luminosities of all flavors are as- 
sumed to be equal, and the RMS average neutrino en- 
ergies are taken to be 12, 22, and 34 MeV, for electron, 
anti-electron, and the other flavors of neutrinos, re- 
spectively. The temperature at the surface of the neu- 
tron star is determined so that heating and cooling by 
neutrino interactions are in equilibrium. The surface 
of the neutron star is arbitrarily defined as the point 
at which density has dropped to p — 10^*^ g cm~'^. 
The equation of state for the electron and positron gas 
includes not only the relativistic pairs as in Paper I 
but also the partially relativistic pairs, which are im- 
portant during the a-process. Once M , R, and Li, are 
specified along with the boundary condition M, we 
obtain numerically the velocity and thermodynamic 
quantities of the neutrino-driven wind as functions of 
r. 

In Paper I, M for each neutrino-driven wind was 
fixed such that the temperature at r = 10, 000 km 
was 0.1 MeV. However in that case no physical solu- 
tion exists for L^, ^ 10^^ ergs s~^, because the tem- 
perature does not cool to 0.1 MeV before the wind 
reaches r = 10, 000 km for any physical solutions. 
Note that equations (l)-(3) have no physical solution 
for M > Mmax- In order to avoid this discontinu- 
ity, we instead adopt M — Mmax, which allows for a 
supersonic wind. The advantage is that the wind so- 
lution exists for any even larger than lO^'^ ergs s~^. 

We evaluate the key parameters for the r-process. 
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Fig. 1. — (a) The maximum mass ejection rates, (b) entropy 
per baryon, and (c) dynamical timescales measured as the time 
for material to cool from T = 0.5 MeV to 0.2 MeV for models A 
(dot-dashed line), B (dashed line), and C (thick-solid line), 
as functions of Ly for the supersonic winds. The thin-solid 
lines in the panels (b) and (c) are for the subsonic wind with 
M = 0.995 X Mmax for model C (see text). Also denoted 
are the results from Sumiyoshi et al. (2000, filled squares and 
circles) and Qian & Woosley (1996, open squares and circles). 
In both cases, the squares and circles are the results with the 
same model parameters M/Mq and R as the models A and B, 
respectively 



i.e., entropy and dynamical timescale, as a function of 
Afmax instead of arbitrary M for each L^. When nu- 
cleosynthesis calculations are performed, we actually 
take M to be slightly smaller than Mmax to obtain so- 
lutions for subsonic winds, as discussed below. As ini- 
tial conditions, we make the reasonable assumptions 
that the electron fraction Yg is set to 0.5, and that no 
Qf-particles or heavy elements have yet formed. 

In this study, we explore three specific models of 
the neutrino-driven winds based upon combinations 
of two parameters M and R. The models given in Ta- 
ble 1 correspond to: A) (M, i?) = {IAMq, 10 km); B) 
(M,i?) = (1.4M0,7km); and C) (M, i?) = (2.0Mq, 10 km). 



The difference between models A and B is the com- 
paction ratio GM/c^R, while neutron star mass is the 
same. On the other hand, effects of general relativity 
are manifest through the compaction ratio. Models 
B and C therefore involve different masses and radii, 
while the compaction ratio is the same. The com- 
paction ratio for models B and C is rather large ~ 0.3. 
However, this value is still consistent with relatively 
soft equations of state for high density matter, or as 
a metastable state pr i or to collapse to a bla ckhole 
( [Brown fc Bethe 1994| ; [Baumgarte et al. 1996|) . Note 
that the case of a smaller compaction ratio <; 0.2 is 
omitted i n this study, where r-p rocessing may not be 
possible ( ICardaU fc Fuller 1997| ). 



Figure 1 shows: a) the mass ejection rate Mmax; 

b) the entropy per baryon S/k &iT ^ 0.5 MeV; and 

c) the dynamical timescale r measured as the time 
for material to cool from T = 0.5 MeV to 0.2 MeV, 
as a function of L^j for supersonic winds. The mass 
ejection rate is the maximum value for physical so- 
lutions of the winds. The dot-dashed, dashed, and 
thick-solid lines respectively denote models A, B, and 
C. As can be seen in Figure 1, for larger L^, Mmax and 
S take larger and smaller values, respectively. A sim- 
ilar trend in entropy can be seen in r. However, the 
time scale saturates at Ljy ~ a few times 10^^ ergs s~^. 
This can be traced to the smaller transonic radii 
and higher initial temperature for higher neutrino lu- 
minosities. The temperature gradient is steeper when 
r is closer to the neutron star surface, while being ap- 
proximately constant for r > (Paper I). Thus, the 
temperature decreases more slowly for r > Tc than for 
r < Tc- For higher luminosity cases J> 10^^ ergs s"""^, 
the winds pass Tc with T > 0.2 MeV. As a result, 
the dynamical timescale does not decrease for higher 
luminosities. 

The r-process favors higher entropy and shorter 
dynamical timescale. Thus, a robust /'-process is dif- 
ficult to obtain once the neutrino luminosity is higher 
than a few times 10^^ ergs s~^. Comparing the re- 
sults of models A and B (or models A and C), we 
find that more "compact" proto-neutron stars eject 
less material, and obtain significantly higher entropies 
and shorter dynamical timescales. This is in good 
agreement with the results of Cardall & Fuller (1997). 
Clearly, the more "compact" neutron star models (B 
and C) will obtain more favorable physical conditions 
for the r-process. 

On the other hand, comparing the results of mod- 
els B and C which are at the same compaction ra- 
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tio, we see that a more massive (or smaller radius) 
proto-neutron star ejects slightly less material and 
provides higher entropy, but involves a larger dynam- 
ical timescale. These can be explained by the rela- 
tions M cx i?5/3^-2^ g ^ R-'^/^M, and r tx RM 
for fixed energy and neutrino luminosity ( |Qian fc 
Woosley 1996) in the Newtonian limit. Based on these 
results alone it is difficult to judge which of the mod- 
els B or C is more favorable for the r-process. Dif- 
ferences of the nucleosynthesis results between these 
models are to be discussed in the next section. 

Figure 1 also shows the results of the numerical 
simulation of the winds from Sumiyoshi et al. (2000) 
for (M, R) = (1.4Mq, 10 km) and = {2.0Mq, 10 km), 
respectively. These are denoted by filled squares 
and circles. Open squares and circles show the re- 
sults of Qian & Woosley (1996) calculated with post- 
Newtonian corrections for the same parameter sets 
of {M,R). Our results for M, S, and r are in 
good agreement with theirs. Ultimately, mass ejec- 
tion rates must be determined by a detailed hydro- 
dynamic simulation of a proto-neutron star with neu- 
trino transport being taken into consideration. Nev- 
ertheless, our mass ejection rates, entropies, and dy- 
namical timescales are in good agreement with those 
of Sumiyoshi et al. (2000), who calculated steady- 
flow neutrino-driven winds with the use of a one- 
dimensional hydrodynamic code including neutrino 
heating and cooling processes. Thus, we conclude 
that the true mass ejection rates are probably close 
to the supersonic winds obtained here. 

Without imposing an outer boundary, however, the 
temperature and density diminish too quickly in the 
supersonic winds. This is not favorable for the r- 
process. Therefore, a slightly smaller mass ejection 
rate, 0.995 x Mmax, is adopted for application to r- 
process calculations in the next section. The advan- 
tage of this reduction is that the temperature and 
density fall off slowly during the r-process (Tg <^ 2.5) 
in these subsonic winds, ensuring enough time for the 
nucleosynthesis of heavy elements to proceed. This 
situation resembles the wind with an outer boundary, 
as if a supersonic wind interacts with a wall shock to 
sustain the temperature and density suitable for the 
r-process. The entropy per baryon and the dynami- 
cal timescale for model C with this M are displayed 
in Figures lb and Ic by thin-solid lines. This small 
reduction of M increases the entropy per baryon only 
about a few percents. The dynamical timescale also 
increases by a factor of < 2 for L^, < 10^^ ergs s~^. 



For Li, > 10^^ ergs s"-'^, however, the dynamical 
timescale increases drastically. This is because the 
winds pass Vc with T > 0.2 MeV for higher luminos- 
ity cases. Furthermore, the velocity and the radial 
gradient in temperature significantly decrease after 
the wind crosses rc for a subsonic wind (Paper I). 
This situation that the dynamical timescale becomes 
significantly longer at an early stage is similar to the 
results of previous detailed hydrodynamic simulations 
([Woosley et al. 1994|; [Witti et al. 1994|). Obviously, 



the r-process would not take place in the winds for 
Li, > 10^^ ergs s~^ owing to the too long dynami- 
cal timescales, as to be confirmed by nucleosynthesis 
calculations in the next section. Note that the en- 
tropy per baryon for Li, > 10^^'^ ergs and the 
dynamical timescale for L^, > 10^^'^ ergs s~^ are not 
displayed in Figures la and lb, since the temperature 
at 10** km, where the calculation is stopped, is higher 
than 0.5 MeV and 0.2 MeV, respectively. 



3. THE r-PROCESS 

Nucleosynthesis yields from the r-process were cal- 
culated for the three neutrino-driven wind models 
discussed in the previous section. We used a fully 
implicit nuclear reaction network for both the a-rich 
freezeout and the subsequent r-process. The network 
consists of over 3000 isotopes all the way from neutron 
and proton up to the plutonium isotopes. We include 
all relevant reactions, i.e., (^,7), (p, 7), (a, 7), {p,n), 
(a,p), {a,n), and their inverse. Reaction rates were 
taken from Thielemann (1995) for nuclei with Z < 46 
and from Cowan, Thielemann, & Truran (1991) for 
those with Z > 47. Weak interactions such as /3- 
decay, /9-delayed neutron emission, electron capture, 
and the capture of electron and anti-electron neutri- 
nos on free nucleons were also included. These neu- 
trino capture processes are of importance as they re- 
duce the number of free neutrons by the "a-effect" 
(Meyer, McLaughhn, & Fuller 1998). Note that fis- 
sion reactions are not included. 

Nucleosynthesis calculations were carried out in 
nine trajectories with neutrino luminosities between 
logL^(ergs s~^) = 52.6 and 51.0 in intervals of 0.2 dex. 
The mass ejection rates were taken to be 0.995 x Mmax 
as described in the previous section. Each calculation 
started when the temperature decreased to Tg = 9 
(where Tg = T/IO^ K). At this point the nuclear 
statistical equilibrium (NSE) consists mostly of free 
nucleons. The initial mass fractions of neutrons and 
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Fig. 2. — Final abundances in model C as a function of atomic mass number, for logL^ (ergs s ^) 
(d), 51.8 (e), 51.6 (f), 51.4 (g), 51.2 (h), 51.0 (i). 



52.6 (a), 52.4 (b), 52.2 (c), 52.0 



protons were therefore given by NSE plus charge neu- 
trality: Xn = 1 — Ye and Xp = Ye, respectively. The 
initial electron fraction Yg was set to be 0.40, which is 
near the lower bound for neutrino-heated ejecta from 
the models of Woosley et al. (1994). This approxi- 
mation is adequate for our study, although ultimately 
the initial electron fraction should be determined by 
the detailed balance of the energy and luminosity be- 
tween electron and anti-electron neutrinos with a hy- 
drodynamic study. 

Figure 2 illustrates the final abundance of each tra- 
jectory for model C as a function of mass number. It 
is found that the trajectories with log (ergs s~^) < 
52.0 (Figure 2d-i) result in a robust r-process. This 
is due to the sufficiently high entropies, S ^ 140fc, 
and short dynamical timescales, r <^ 30 ms (Fig- 
ure 1). Figure 3 shows the neutron-to-seed ratio 
Yn /Ysccd (solid line), the mass fraction of seed nuclei 
Xseed (dashed line), and the electron fraction Ye (dot- 
dashed line) at Tg = 2.5, which is approximately the 
temperature at the beginning of the r-process phase. 



is also shown by the dotted line. As can be seen in 
Figure 3, the neutron-to-seed ratio reaches the max- 
imum ~ 160 at logL^(ergs s~^) = 51.8. As a result, 
for the corresponding trajectory, the third abundance 
peak and the heavier are prominent (Figure 2e) . How- 
ever, the total abundance of r-process nuclei produced 
from this trajectory is small. This is because the mass 
fraction of seed nuclei obtains the minimum value at 
this luminosity, as can be seen in Figure 3. Note 
that Xr is smaller for > 10^^ ergs s~^ despite 
higher Xgccd owing to the smaller neutron-to-seed ra- 
tio. For lower luminosities log Li^(ergs s""'^) < 51.6, 
/9-delayed neutron emission smoothes the abundance 
pattern as seen in Figure 2, because of the longer dy- 
namical timescales (Figure 1). 

As anticipated in the previous section, trajectories 
for logLjy > 52.2 (Figures 2a-c) do not contribute sig- 
nificantly to the r-process, because of the too long dy- 
namical timescales (Figure Ic). Neutrino capture on 
free nucleons somewhat increases the electron fraction 
to Ye ~ 0.42 - 0.43 at Tg = 2.5 by the a-effect (|Meyer| 



Heie, "seed" refers to all nuclei with A> 12. The fi- 



nal mass fraction of r-process elements {A > 100) Xr 



et al. 1998) as seen in Figure 3 (dot-dashed line). 



Thus, it is important to include the neutrino capture 
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Fig. 3. — The neutron-to-seed ratio (solid line), the mass 
fra<;tion of seed nuclei (dashed line), and the electron fra<;tion 
(dot-dashed line) at Tg = 2.5. The final mass £ra<;tion of r- 
process elements (A > 100) is also denoted by the dotted line. 



reactions in order to correctly calculate the resulting 
r-process nucleosynthesis. This effect is not, however, 
as serious as that suggested by Meyer et al. (1998) in 
the present r-process conditions. The reason is that 
the dynamical timescales in our models r 30 ms 
are sufficiently shorter than the neutrino interaction 
lifetime for the neutrino luminosities considered. The 
dynamical timescale of Meyer et al. (1998) was fixed 
at a value as large as t = 0.3 s. 

In order to compare the nucleosynthesis results 
of each model with the solar-system r-process abun- 
dances, the yields were integrated over the mass- 
weighted time history in the following manner. The 
time evolution of the neutrino luminosity is approxi- 
mately given by 



= i,.oexp(-i/ry), 



(5) 



where ij^o is the initial neutrino luminosity at the 
neutrino sphere, and is the cooling timescale. We 
take L^Q to be w 4 x 10^^ ergs s~^ (logL^o = 52.6) 
from Wooslcy et al. (1994). In principle L^q depends 
on the radius of neutrino sphere and should be cor- 
rected for the gravitational redshift effect on neutrino 
energy. However, the present approximation is ade- 
quate for our purposes as neutrino-driven winds with 
4 X 10^^ ergs s~^ do not substantially contribute 
to the r-process yields. The cooling timescale is set 
equal to 1.25 s so that the total time integrated neu- 
trino energy becomes 3 x 10^^ ergs. With this t^, the 
Figures 2a-i correspond to time slices at 0, 0.6, 1.2, 
1.7, 2.3, 2.9, 3.5, 4.0, and 4.6 s, respectively. 

We integrate the mass-weighted r-process yields 
assuming that the trajectories can be described at 
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Fig. 4. — The mass-weighted integrated yields for models A 
(a), B (b), and (c) as functions of mass number (lines). Also 
denoted are the scaled solar r-abundances (points). 



each time by steady flow corresponding to the neu- 
trino luminosity and neutron-star radius at that time. 



where 



M(t)dt 



(6) 
(7) 



is the total mass ejected from the neutron star. 

Obviously, our assumption of steady flow for each 
trajectory is oversimplified. However, the thermal 
properties of material in the neutrino-heated wind are 
mostly determined near the surface of the neutron 
star ^ 30 km, where the equilibration time 0.1 s) 
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is shorter than the neutrino luminosity timescale 

1 s). In addition, during the relaxation of the 
neutron star radius ~ 2 s, the neutrino luminosity 
is so high Li, ^ 10^^ ergs s~^ that few r-process el- 
ements are synthesized. Hence, the assumptions of 
steady flow for each neutrino luminosity and a fixed 
neutron-star radius are probably reasonable. 

In Figure 4, the integrated yields for models A 
(a), B (b), and C (c) are compared with the solar- 
system r-process abundances (Kappeler, Beer, & Wis- 
shak 1989). In models B and C the solar r-process 
abundances are normalized to the third peak. For 
model A the second peak is used for normalization. 
Comparison of the three results suggests that the 
"compaction" of the proto-neutron star is essential 
to reproduce the pattern of the solar r-process abun- 
dances. Indeed, the abundance pattern for elements 
with A > 120 in models B and C are in excellent 
agreement with the scaled solar r-abundances, while 
only a weak r-process occurred in model A. 

In all three models there is about a factor of 10 
overproduction of nuclei with A <^ 120. For A 90 
in particular, there is a prominent peak as also seen 
in previous studies ( Freiburghaus et al. 1999). In 
the present models, however, the overproduction is 
about a factor of 10 smaller than that of Woosley et al. 
(1994). Moreover, if the abundances of A - 100-130 
were smoothed, they would agree reasonably well with 
the solar r-pattern. 

The good agreement in models B and C is a con- 
sequence of the robust r-process which starts at the 
early phase 2 s) of the neutrino-driven wind. At 
this time is still as high as ^ 10^^ ergs s~^ and the 
integrated yields are dominated by the matter ejected 
there. This is in contrast to the result of Woosley et 
al. (1994) in which the r-process occurs only during 
the late phase of the neutrino-driven wind 10 s), 
where the mass ejection rate has already significantly 
declined. 

In Table 2, the total ejected mass Mcj and mass of 
r-process {A > 100) elements Afcj, r for each model are 
compared to those of Woosley et al. (1994). Note that 
Mej^r in the W94 column is the sum of the last 16 tra- 
jectories in Woosley et al. (1994), whose pattern was 
in excellent agreement with the solar r-abundances. 
The ejected r-process mass in model C is about a fac- 
tor of two larger than that of Woosley et al. (1994) 
owing to production of r-nuclei in early times when M 
is high. It is interesting to note that the total ejected 
masses of all models are smaller than that of Woosley 




100 



150 
mass number 



200 



Fig. 5.- 

{1.7Af0,lO 



— Same as the Figure 4, but for (M, R) 
km), (1.8Mq, 10 km), and {1.9Mq, 10 km). 



et al. (1994), while M is set to be nearly equal to the 
maximum value in this study. This is a consequence of 
our more compact neutron stars, which obtain smaller 
M as can be seen in Figure la. In particular, the ra- 
dius of the neutron star in Woosley et al. (1994) is 
much larger than 10 km in the first a few seconds, 
while those in this study are taken to be constant. In 
addition, M in Woosley et al. (1994) during the first 
~ 1 s is significantly larger than expected for a steady 
state wind. 

Comparing the three models in Table 2 we see 
that Mcj does not significantly change from model 
to model, being somewhat smaller in models with 
a deeper gravitational potential. As can be seen in 
Figure 4, however, the r-abundance pattern is sub- 
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Table 2 
Total Ejected Mass 



A 


B 


C 


W94t 


M^-JMq 2.9 X 10-3 


1.7 X 10-3 


1.1 X 10-3 


3.2 X 10-2 


Mej,r/M0 5.3x10-* 


2.0 X 10-5 


1.4 X 10-* 


5.8 X 10-5* 


tWoosley et ah (1994) 



t abundances produced in the last 16 trajectories 



stantially different among these models. Model A re- 
produces the first {A ~ 80) and second {A ~ 130) 
r-proccss peaks, but no elements with A 130 are 
synthesized. In contrast, the solar r-pattern, includ- 
ing all three r-process peaks, is well reproduced in 
the "compact" neutron star models B and C. This 
is a consequence of the significantly higher entropies 
S ^ ISONAk and shorter dynamical timescales r ^ 
30 ms for these models. 

Comparing the absolute yields in models B and 
C, we find that the neutrino- driven wind of a more 
massive proto-neutron star produces more r-process 
elements when the compaction ratio is the same. In- 
deed, the mass fraction of r-elements in model C is 
about a factor of 10 larger than that of model B. The 
reason is that the dynamical timescale in model B 
is shorter than in model C (cf. Figure 1). In such 
fast winds, fewer seed nuclei are produced, although 
the neutron-to-seed ratio is high enough for a robust 
r-process (see also Hoffman, Woosley, & Qian 1997). 

Let us discuss how sensitive these results are to 
the compaction ratio. We present in Figure 5 the 
mass-weighted integrated yields for three more cases 
(M,i?) = (1.7M0,1O km), (1.8M0,1O km), and 
(1.9M0, 10 km). These correspond to compaction ra- 
tios of 0.25, 0.27, and 0.28, respectively. Interest- 
ingly, the pattern in the first model is not signifi- 
cantly different from that of model A. In fact, it is 
in slightly better agreement with the first and second 
solar r-process peaks but does not produce heavier 
elements (Figure 5a). The second model produces a 
small amount oi A > 130 elements, and the third 
one reasonably reproduces the solar r-pattern up to 
the third peak (Figures 5b and c) . These calculations 
indicate that the r-process strongly depends on the 
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Fig. 6. — Comparison of the mass radius relation (lines) for 
various equations of state with the models (circles) A, B, and 
C considered in this work. See text for an explanation of the 
various equations of state. 

compaction ratio. A good fit is only obtained when 
the compaction ratio is ^ 0.28. 

4. DISCUSSIONS AND CONCLUSION 

We have explored the r-process in neutrino-heated 
ejecta (neutrino-driven winds) from nascent neutron 
stars for various neutron-star masses, radii, and neu- 
trino luminosities. Our results suggest that proto- 
neutron stars with a large compaction ratio (GM/c^R ~ 
0.3; models B and C), where general relativistic ef- 
fects are important, provide the best conditions for 
the r-process. Such compact proto-neutron stars lead 
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to sufficiently high entropy S ^ 100 — 200NAk and 
short dynamical timescales t <^ 30 ms to produce 
an optimum neutron to seed ratio for material in 
the neutrino-driven winds. The total time-integrated 
yields from these models are in good agreement with 



neutrino emission or a late softening and heating by 
kaon condensation (Thorsson, Prakash, & Lattimer 
1994). For example, the softening of a proto-neutron 
star by kaon condensation is estimated to occur in a 



the solar r-process abundance pattern especially for 



nuclei with A ~ 120 — 200. Our short dynamical 
timescale and assumed condition of steady-flow equi- 
librium for material in the winds lead to a good r- 
process even at early times when the neutrino lumi- 
nosity is still high L^, ~ 10^^ ergs s~^. As a result, 
the overproduction of nuclei A ^ 120 is significantly 
reduced, compared to that of Woosley et al. (1994), 
although being still overproduced by a factor of ~ 10. 

Our results also suggest that, for models of the 
same compaction ratio, e.g., B and C, a higher neutron- 
star mass yields more r-process material. This is due 
to the somewhat longer dynamical timescale for ma- 
terial ejected from high-mass neutron stars which in- 
creases the amount of seed nuclei. 

The implication of these two results is that an op- 
timum r-process which does not overproduce A 90 
nuclei may require both a well developed bubble and 
a massive proto-neutron star while the r-process el- 
ements are synthesized. It should be noted that we 
presume that the bubble is already evacuated and a 
steady-state wind developed for all neutrino luminosi- 
ties. For the dynamical supernova model of Woosley 
et al. (1994), however, the proto-neutron star is still 
relaxing and the bubble slowly expanding during the 
first ^ 2 s. In the model C of this study, the r- 
process starts at — 10^^ ergs (Figure 2), which cor- 
responds to 1.7 s in equation (5). At L,j ~ lO'^^ ergs, 
therefore, the entropy per baryon might be somewhat 
overestimated, although the dynamical timescale as 
a function of Liy shows a similar trend to those of 
previous hydrodynamic studies (Figure Ic). A de- 
tailed hydrodynamic study of the neutrino-induced 
explosion for a massive supernova progenitor will be 
eventually needed to confirm whether the early stage 
of the neutrino-driven winds from compact neutron 
stars are really the astrophysical site of the r-process. 

Alternatively, this scenario might occur if the low 
entropy material initially ejected from the supernova 
were to fall back onto the neutron star. This would 
produce the massive remnant. This epoch might then 
follow by a second epoch of high neutrino flux into the 
previously evacuated bubble. This might occur, for 
example, if the proto-neutron star were to experience 
a late transition to strange matter and subsequent 



few to 10 seconds after core bounce (Keil & Janka 
T995| ), which may be relevant to the r-process in this 



study. In either case, it is clear that these new results 
could suggest an interesting new twist to the super- 
nova paradigm. 

If these models for the r-process are correct, then 
our results could also pose a significant constraint on 
the equation of state for the neutron-star remnant. 
A summary of the mass-radius relationship of some 
contemporary equations of state is shown in Figure 6. 
These are compared with the various models consid- 
ered in this work. These representative equations 
of state include an example of a non-relativistic nu- 
cleon potential model equation of state with both two- 
body and three-body terms (AV14 + UVII; Wiringa, 
Fiks, & Fabrocini 1988) as well as various relativistic 
mean field equations of state, including those based 
upon a relativistic Hartree (G; Glendenning 1989 and 
HV; Glendenning 1985), Hartree Fock (HVF; Weber 
& Weigel 1989) calculation and those with strange- 
matter interiors (SSI & SS2; Glendenning & Weber 
1992). 

Model A with M = IAMq and i? = 10 km is a 
typical neutron star mass and a typical radius for a 
non-relativistic EOS and several relativistic ones as 
well. However, it does not produce a good r-process 
in our simulations. Generating the observed abun- 
dance curve seems to require a slightly more exotic 
paradigm. If we consider our best r-proccss model, 
i.e. model C with M = 2.0Mq and i? 10 km, 
we see that many of the equations of state lead to 
such a remnant. In all cases, however, such a mas- 
sive remnant is very close to or equal to the maxi- 
mum mass allowed by the EOS. If such a remnant 
is required to produce the r-process, then it seems 
that the r-process must occur for a remnant that has 
accreted to its maximum mass and is about to col- 
lapse to a black hole. This suggests another possible 
paradigm in which the r-process occurs in the pre- 
viously evacuated bubble by neutrinos generated as 
the short-lived proto-neutron star begins its collapse 
to a black hole (see also Qian, Vogel, & Wasserburg 
1998). This would perhaps require that the r-process 
occurs in a more massive progenitor star which does 
not leave a neutron star remnant such as probably 
occurred in SN 1987A (|Brown fc Bethc 1994]). 
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Although model B with M = IAMq and i? = 7 km 
also gives a reasonable r-process abundance distribu- 
tion, it is very difficult to form such a compact neu- 
tron star of this mass. Such lower-mass high com- 
paction stars could perhaps form if a strange matter 
core develops as in SSI and SS2 of Figure 6. 

The results of this study are also of importance 
for the chemical evolution of the Galaxy. Ishimaru & 
Wanajo (1999) have suggested that the large disper- 
sion nf eurnpiiuTi with respect to iron amnnp; metal- 



poc 



r sta,rs in the Cala.ctic halo ( |Rya,n et a,1. 199f| : 



McWilliani 1997) can be reproduced if the r-process 
elements originate from stars of > SOA/q or 8 — IOMq. 
Our results also imply that, for a given equation of 
state, a more massive proto-neutron star leads to a 
smaller compaction ratio and therefore more favor- 
able conditions for the r-process. Furthermore, if the 
maximum mass of a proto-neutron star were as large 
as ^ 2 A/0, the progenitor star is probably very mas- 
sive. Such stars with > 30Mq account for only ^ 10% 
of all supernova progenitors with a typical initial mass 
function. Nevertheless, the amount of r-process ma.t- 
ter in model C is about a factor of 2 higher than that 



of Woosley et al. (1994), so that such a restriction to 



high mass stars may be reasonable. 

It should be noted that, even if neutrino-driven 
winds fail to reproduce the solar r-pattern espe- 
cially of the third peak, they still could be signifi- 
cant sources of the solar r-elements of at least the 
second peak and lighter. If that is the case, heav- 
ier r-process elements must be synthesized in other 
sources. This is consistent with the scenario of 
mu ltiple r-process sites implied bv meteoritic abun- 



da nces fWasserburg. Busso. fc Gallino 1996). spec- 



troscopic studies of metal-poor halo stars (Snedcn et 
al. I 2000), and the chemical evolution studies (Ishi- 
Wanajo 200C1| ). This could be also explained 



by the neutrino-driven wind models described here 
which show the strong dependence of the r-process 
yields on the remnant mass. 

In fact, most of the metal-poor stars in the Galac- 
tic halo show the higher ratios of light-to-heavy r- 
elements (e.g., Sr/Ba or Sr/Eu) than that of the solar 
r-abundances, although the values differ f rom star to 
star b y up to 2 orders of ma gnitude ( Mc William 
19c(|; [Ishimaru fc Wanajo 2000[) . This implies that 
these r-elements might originate from the "weak" r- 
processing associated with smaller compaction ratios 
< 0.28 (Figures 4 and 5). However, some metal-poor 
stars (e.g., CS 22892-052) show somewhat smaller ra- 



tios of light-to-heavy r-elements compared to that of 
the solar r-abundances, by a factor of ~ 2 — 3 (Sneden 
et al. 2000; see also Wasserburg et al. 1996 for a me- 
teoritic study), which cannot be reproduced by any 
models in this study. If such low light-to-heavy ratios 
originate from neutrino-driven winds, we need some- 
what more compact proto-neutron stars than consid- 
ered in this study, otherwise we have to seek other 
causes to reduce the lighter r-nuclei, e.g., fallback of 



the fi rst, low entropy material ejected (Woosley et al 



19941 ) 



Obviously, our approach of assuming constant val- 
ues for the initial electron fraction, neutrino lumi- 
nosity in each trajectory, the radius of the proto- 
neutron star, as well as an exponential time evolution 
of neutrino luminosity, is too simplified to conclude 
that these neutrino-driven wind models are an accu- 
rate description of the r-process. For example, the 
time variation of Ye is of particular importance in 
the r-process. Recent hydrodynamic studies of core- 
collapse supernovae with accurate treatments of the 
neutrino transport have shown that Yg exceeds 0.5 for 
the first a few hundred ms after core bounce ( Rampp 
fc Janka 2GG0| ). If that is the case, then the overpro- 
duction of A « 90 nuclei in our results would disap- 



pear (Hoffman ct al. 1996). In the present study, the 
initial electron fraction was fixed at a relatively low 
value, 0.4, although it increases to 0.42 — 0.43 by the 
time of onset of the r-process. However, Ye may be 
significantly higher than that during the early phase 
of the neutrino-driven wind. Even with a higher Ye, 
however, a sufficient neutron-to-seed ratio could be 



obtained if the a-process is less efficient ( Hoffman et 
al. 1997). Instead, the yields of r-process elements 



might decrease owing to a smaller amount of seed 
material. 

Future hydrodynamic studies of core-collapse su- 
pernovae with accurate neutrino transport, including 
multidimensional effects such as convection and ro- 
tation, as well as magnetic fields will probably ulti- 
mately be required before we can be certain that the 
neutrino-driven winds are indeed the major r-process 
site in the Galaxy. Nevertheless, our results have con- 
firmed that the scenario of neutrino-driven winds is 
still viable and promising as the true astrophysical 
site for the r-process. 

We would like to acknowledge useful discussions 
with Y. Ishimaru, N. Itoh, K. Sumiyoshi, and M. 
Terasawa. This work was supported in part by Japan 



11 



Society for Promotion of Science, and by the Grant- 
in Aid for Scientific Research (1064236, 10044103, 
11127220, and 12047233) of the Ministry of Educa- 
tion, Science, Sports and Culture of Japan, and DoE 
Nuclear Theory grant DE-FG02-95-ER40394. 

REFERENCES 

Baumgartc, T. W., Janka, H. -T., Koil, W., Shapiro, S. 
L., & Tcukolsky, S. A. 1996, ApJ, 468, 823 

Brown, G. W. & Botlic, H. A. 1994, ApJ, 423, 659 

Cardall, C. Y. & Fuller, G. M. 1997, ApJ, 486, Llll 

Cowan, J., Thielemann, F. -K., & Truran, J. W. 1991, 
Phys. Rep., 208, 267 

Preiburghaus, C., Rembges, J. -P., Rauscher, T., Kolbe, 
P., Thielemann, P. -K., Kratz, K. -L., Pfeiffer, B., & 

Cowan, J. J. 1999, ApJ, 516, 381 

Preiburghaus, C, Rosswog, S., & Thielemann, P. -K. 
1999, ApJ, 525, L121 

Glendenning, N. K. 1985, ApJ, 293, 470 

Glendenning, N. K. 1989, Nucl. Phys. A, 493, 521 

Glendenning, N. K. & Weber, P. 1992, ApJ, 400, 647 

Hoffman, R. D., Woosley, S. E., Fuller, G. M., & Meyer, 
B. S. 1996, ApJ, 460, 478 

Hoffman, R. D., Woosley, S. P., & Qian, Y. -Z. 1997, ApJ, 
482, 951 

Ishimaru, Y. & Wanajo, S. 1999, ApJ, 511, L33 

Ishimaru, Y. & Wanajo, S. 2000, in Ftrst Stars, A. Weiss, 
T. Abel,&: V. Hill eds. Springer- Vorlag, Berlin, p. 189 

Itoh, N., Hayashi, H., Nishikawa, A., & Kohyama, Y. 
1996, ApJS, 102, 411 

Kappeler, P., Beer, H., & Wisshak, K. 1989, Rep. Prog. 

Phys., 52, 945 

Keil, W. & Janka, H. -Th 1995, A&A, 296, 145 
McWilliam, A. 1997, ARA&A, 35, 503 
McWilliam, A. 1998, AJ, 115, 1640 

Meyer, B. S., Mathews, G. J., Howard, W. M., Woosley, 
S. E., & Hoffman, R. D. 1992, ApJ, 399, 656 

Meyer, B. S., McLaughlin, G. C, & Fuller G. M. 1998, 
Phys. Rev. C, 58, 3696 



Otsuki, K., Tagoshi, H., Kajino, T., & Wanajo, S. 2000, 

ApJ, 533, 424 (Paper I) 

Qian, Y. -Z. & Woosley, S. E. 1996, ApJ, 471, 331 

Qian, Y. -Z., Vogel, P., & Wasserburg, G. J. 1998, ApJ, 
494, 285 

Qian, Y. -Z. 2000, ApJ, 534, L67 

Rampp, M. & Janka, H. -T. 2000, ApJ, 539, L33 

Ryan, S. G., Norris, J. P., & Beers, T. C. 1996, ApJ, 471, 
254 

Shapiro, P., & Tcukolsky, S. 1983, Black Holes, White 
Dwarf's, aud Neutron Stars (New York: Wiley) 

Sneden, C, McWilham, A., Preston, G. W., Cowan, J. J., 
Burris, D. L., & Armosky, B. J. 1996, ApJ, 467, 819 

Sneden, C, Cowan, J. J., Debra, L. B., & Truran, J. W. 

1998, ApJ, 496, 235 

Sneden, C, Cowan, J. J., Ivans, 1. 1., Puller, G. M., Buries, 
S., Beers, T. C., Lawler, J. E. 2000, ApJ, 533, L139 

Sumiyoshi, K., Suzuki, H., Otsuki, K., Terasawa, M., &: 
Yamada, S. 2000, PASJ, 52, 601 

Takahashi, K., Witti, J., & Janka, H. -Th. 1994, A&A, 

286, 857 

Thielemann, F.-K. 1995, private communication 

Thorsson, V., Prakash, M., & Lattimer, J. M. 1994, 
Nucl. Phys. A, 572, 693 

Tsujimoto, T., Shigeyama, T., & Yoshii, Y. 2000, ApJ, 
531, L33 

Wasserburg, G. J., Busso, M., & Gallino, R. 1996, ApJ, 
466, L109 

Weber, P. & Weigel, M. K. 1989, Nucl. Phys. A, 505, 779 

Wiringa, R. B., Piks, U., & Pabrocini, A. 1988, 
Phys. Rev. C, 38, 1010 

Wheeler, J. C, Cowan, J. J., & Hillebrandt, W. 1998, 

ApJ, 493, LlOl 

Witti, J., Janka, H. -Th., & Takahashi, K. 1994, A&A, 
286, 841 

Woosley, S. E. & Hoffman, R. D. 1992, ApJ, 395, 202 

Woosley, S. P., Wilson, J. R., Mathews, G. J., Hoffman, 
R. D., & Meyer, B. S. 1994, ApJ, 433, 229 



This 2-column preprint was prepared with the AAS lAT^^X 
macros v4.0. 



12 



